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ABSTRACT 

It  is  shown  that  the  interaction  between  macroscopic,  non- 
dissipative  media,  and  time-varying  electromagnetic  fields  can  be 
described  by  a  time-averaged  potential  function.  From  this  function 
it  is  possible  to  derive  phenomenologically  the  tensors  that  describe 
any  of  the  us  vial  electro-  and  magneto-optic  effects  for  electric  and 
magnetic  fields  of  any  frequency.  In  addition,  these  same  potential 
functions  describe  the  various  optical  nonlinearities  like  harmonic 
generation  in  KDP,  and  harmonic  generation  by  electric  quadrupole 
and  magnetic  dipole  nonlinearities. 

The  symmetry  relations  first  derived  by  Armstrong,  Bloembergen, 
Ducuing,  and  Pershan  for  electric  dipole  nonlinearities  follow  directly 
from  the  methods  presented  here.  In  addition,  one  can  derive 
analogous  relations  for  electric  quadrupole  and  magnetic  dipole 
nonlinearities.  These  relations  also  demonstrate  the  reciprocal 
nature  of  the  linear  electro-optic  effect  and  rectification  of  light. 

The  Faraday  effect  and  the  production  of  a  dc  magnetization  due  to 
incident  circularly  polarized  light  are  also  reciprocal  effects. 
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Nonlinear  Optical  Properties  of  Solids:  Energy  Considerations 

by 

P.  S.  Pershan 

Division  of  Engineering  and  Applied  Physics 
Harvard  University,  Cambridge,  Massachusetts 

I.  Introduction 

The  nonlinear  responses  of  macroscopic  media  to  applied  electric  and 
magnetic  fields  at  lower  than  optical  frequencies  are  well  known.  Crystal 
rectifiers  and  harmonic  generators,  magnetic  amplifiers,  saturable  reactors, 
etc.  ,  are  just  a  few  of  many  examples. 

Until  recently,  the  only  available  electromagnetic  [  E.  and  M.  ]  fields 
at  optical  frequencies  were  so  feeble  that  one  could  linearize  the  response  of 
any  material  body  to  them  wtih  no  appreciable  error.  With  the  development 
of  optical  lasers,  this  is  no  longer  true  and  recent  experiments  [1,  2,  3,  4, 

5,  6]  have  clearly  demonstrated  nonlinear  responses  to  optical  frequency  fields. 

Theoretically,  these  optical  phenomena  have  been  considered  from  both 
a  quantum  mechanical,  or  first  principle,  approach  as  well  as  phenomenolo- 
gically[7,  8,  9,  10,  11,  12,  13]. 

It  is  the  purpose  of  this  paper  to  show  that  energy  considerations  first 
suggested  by  Armstrong,  Bloembergen,  Ducuing,  and  Pershan  [  ABDP  ]  can 
be  generalized  to  include  all  nonlinear  electromagnetic  processes  in  which 
the  medium  is  nonabsorptive .  In  particular,  in  addition  to  the  ones  first  ob¬ 
tained  by  ABDP  for  electric  dipole  nonlinearities,  it  is  possible  to  derive 
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dispersion  symmetry  relations,  for  electric  quadrupole  and  magnetic  dipole 
optical  nonlinearities,  solely  from  macroscopic  arguments.  The  method  to 
be  presented  is  so  general  that  it  will  be  possible  to  use  one  form  of  tensor 
to  describe  all  nonlinear  effects  of  the  same  type  regardless  of  the  frequencies 
involved.  For  example.,  the  dc  magnetic  Voigt  effect,  microwave  modulation 
of  light  by  magnetic  techniques,  and  magnetic  optical  nonlinearities  are  all 
described  by  different  Fourier  components  of  the  same  tensor.  With  the  aid 
of  these  dispersion  symmetry  relations,  it  will  be  possible  to  estimate  the 
order  of  magnitude  of  some  of  the,  as  yet,  unobserved  optical  nonlinearities. 

The  thermodynamic,  or  energy,  arguments  will  be  developed  in 
Section  II.  Sections  III  through  V  will  be  concerned  with  the  various  specific 
types  of  nonlinearities. 
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II.  Maxwell's  Equations  —  Energy  Considerations 
The  starting  point  from  which  one  derives  the  E.  and  M.  wave  equations 
for  macroscopic  media  is  Maxwell's  equations  in  vacuum  and  a  distribution  of 
charges  and  currents  [14,  15]. 

V  X  e  =  -  (1/c)  ab/at  (2-1) 

V  X  b  =  (l/c)  ae/at  +  (4'Tr/c)  j  . 

Taking  a  suitably  defined  average,  one  obtains 

V  X  g  =  -  (1/c)  a  B/  at  (2-2) 

V  X  B  =  (1/c)  aE/at  +  (4ir/c)  J 

where  [16] 

j  =  ap/at  +  c  V  X  M  -  a(v  •  Q)/at  +  •  •  •  (2-3) 

P  is  the  electric  dipole  moment  per  unit  volume,  M  is  the  magnetic  dipole 
moment  per  unit  volume,  and  Q  is  the  electric  quadrupole  moment  per  unit 
volume.  It  is  well  known  that  the  dipole  moment  per  unit  volume  is  uniquely 
defined  only  when  the  net  charge  density  per  unit  volume  is  zero.  Similarly, 

M  and  Q  are  not  uniquely  defined  when  a  P/  at  and  P  respectively  are  non¬ 
zero.  It  is  possible,  however,  to  define  M  and  Q  in  a  meaningful  manner. 

In  the  usual  manner  one  obtains  the  following  form  of  the  energy  con¬ 
servation  equation. 

(c/4ir)  V  •  (E  X  B)  +  (1/4tt)  B  .  a  B/ a  t 
+  (1/4it)  E  •  aE/at  +  E  •  J  =  0  .  (2-4) 

For  a  vacuum,  i.  e.  ,  J  =  0,  the  first  term  is  the  power  flow,  and  the  second 
and  third  terms  are  the  time  derivative  of  the  energy  density  per  unit  volume. 
The  interaction  between  material  and  E.  and  M.  field  is  thus  represented  by 
the  last  term.  It  is  important  to  realize  that  the  last  term  is  not  simply  the 
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contxibution  of  the  material  to  the  energy  density  per  unit  volume.  Only  when 
M  and  Q  vanish  can  we  regard  E  •  J  =  E  •  d  P/  9t  as  the  time  derivative  of  an 
energy  density.  This  was  the  case,  for  example,  in  ABDP  where  only  the 
electric  dipole  nonlinearity  was  treated.  In  general,  one  must  do  a  partial 
integration  to  obtain 

(c/4it)  V  .  (E  X  H  -  4Trc"^E,-  dQ/dt) 

+  (l/4Tr)H-  8B/at  +  (l/4ir)  E  •  dD/dt 

+  V  E  :  dQ/dt  =  0  (2-5) 

where  H  =  B  -  4Tr  M  and  P  =  E  +  4Tr  p.  Note  that  the  above  definition  of  D 
implies 

V  •  p  +  4Tr  V  V:  5  =  4irp. 

For  the  purposes  of  this  paper,  it  is  not  important  whether  this  definition 
or  one  in  which  V.  D  =  4Trp  <is  given.  The  terms  -  cV  •  (£  iand- 

-V  (E'8^^t)  represent  divergences  of  an  energy  flow  through  the  material 
medium.  This  is  demonstrated  in  Fig.  I  for  the  E  x  M  term.  The  circular 
current  corresponds  to  a  magnetization  out  of  the  page.  The  charge  gains 
energy  from  the  E  field  at  and  gives  it  to  the  E  field  at  z^;  there  is  a  net 
transfer  of  energy  from  left  to  right,  opposite  to  the  vector  cross  product 
E  X  M.  Similar  considerations  will  demonstrate  the  power  flow  for  the 
quadrupole  term. 

If  the  material  energy  density  per  unit  volume  is  U, 

8U/8t=H-  8M/8t  +  E-  8  P/ 8 1  +  V  E;  (8  Q/ 8 1)  (2-6) 

so  that  Eq.  2-5  becomes 

V  •  S  +  (l/4ir)  H  •  8H/8t  +  (l/4ir)  E  •  8E/8t 

+  8U/8t  =  0  (2-7) 

where  S  is  Poynting's  vector  generalized  to  include  quadrupole  effects. 


M  -OUT  OF  PAGE 

O 


Z 


E  X  M 


FIG.  1.  SCHEMATIC  DEMONSTRATION  OF  POWER  FLOW  PROPORTIONAL 
TO  -  (ExM).  THE  COUNTERCLOCKWISE  CURRENT  J_  IS  EQUIVA¬ 
LENT  TO  A  MAGNETIZATION  OUT  OF  THE  PAGE. 
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Equation  2-7  says  that  the  rate  at  which  energy  flows  out  of  a  volume  element 
is  equal  to  the  rate  at  which  the  energy  stored  in  the  E.  and  M.  field  is  de¬ 
creasing  (i.  e.  ,  -  (l/ 4ir)  [H  •  3  H/  3t  +  E  •  9E/  9  t])  plus  the  rate  at  which  the 
material  is  doing  work  on  the  E.  and  M.  field  (i.  e.  ,  -3U/9t).  Equation  2-6 
can  be  recognized  as  the  usual  expression  for  the  work  done  on  a  system  by 
the  external  fields  (17], 

Restricting  our  attention  to  nondissipative  media,  it  is  clear  that  in 
the  steady  state  the  average  work  done  on  (or  by)  the  material  must  be  zero. 

In  the  linear  problem,  for  example,  E  and  3  P/3  t  are  90  degrees  out  of  phase 
and  the  average  of  their  products  vanishes.  This  is  not  to  say,  however,  that 
the  average  energy  stored  in  the  medium  by  virtue  of  its  polarization  is  zero. 
When  the  fields  are  initially  turned  on,  work  is  done  on  (or  by)  the  material 
to  establish  what  is  eventually  termed  the  "steady  state"  amplitudes  of  P,  M, 
Q,  etc.  It  is  most  reasonable  to  assume  that  the  net  work  done  in  producing 
the  steady  state  amplitudes  is  independent  of  the  manner  in  which  they  were 
established.  Equivalently, the  final  "steady  state"  energy  density  only  depends 
on  the  "steady  state"  fields  and  polarizations,  not  on  how  they  were  produced. 

In  order  to  make  these  statements  more  quantitative,  consider 

t  +  T/2 

jg^w,  t)  =  ^  y  ^t)  exp  (-i«t)  dt  (2-8) 

t  -  T/2 

where  T  is  a  time  long  enough  so  that  wT  >>  1,  but  still  short  enough  that 
6(u,  t)  is  independent  of  T.  If  we  were  only  concerned  with  linear  problems, 
this  last  point  would  mean  that  T  is  much  less  than  the  time  in  which  the 
fields  go  from  zero  to  their  "steady  state"  values.  Actually,  we  are  primarily 
interested  in  nonlinear  problems  so  that  T  must  be  small  compared  to  the 
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time  in  which  the  nonlinear  interactions  will  significantly  change  the  Fourier 
distribution  of  any  of  the  field  variables.  This  means  that  we  are  restricted 
to  weak  nonlinearities.  Fortunately,  this  is  not  a  serious  restriction  since 
the  nonlinearities  in  which  we  will  be  interested  are  weak  enough  that  there 
is  no  difficulty  in  picking  a  time  T.  The  Fourier  distribution  of  all  variables 
is  defined  similarly  to  £q.  2-8  and  will  be  indicated  by  script  capitals.  If  there 
are  n  interacting  waves  ;  ^ 

E(t)  =  2  Re  ^  6  ,  t )  exp(iio^t)  (2-9) 

V  =  1 

and  Eq.  2-6  becomes 


n 

8U/8t=2Re  ^  [  3C*  .  t)  •  810  (w^,  t)/ 8  t 

V  =  1 


‘"v  ® 


n 

+  2  Re  Y  [i*  ("v»  t)  •  9£ 
V  =  1 

+  iw^,  J.*  ,  t)  •  iP  (w^,  t)] 


n 

+  2  Re  ^  [V  §♦  (u^,  t)  :  8£(w^,  t)/8t 
V  =  1 

+  i«^  V  t)  («^.  t)] 

+  high  frequency  terms .  (2-10) 

The  physical  statements  we  have  been  able  to  make  about  the  energy  only  per¬ 
tained  to  the  time  average  values.  Averaging  over  T  causes  the  high  frequency 
terms  in  Eq.  2-10  to  vanish. 

Under  steady  state  conditions  STTl  («  /  t)/ 8t  =  0,  9jE  (w  ,  t)/ 8t  =  0, 

dfi  (w.  ,  t)/ 8t  =  0  and  since  one  requires  that  <8  U/dt^,r>  =  ^ 
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n 

0  =  2  Re  I  (u)^.  t) .  n  (cj^,  t)  +  i*  (w^.  t)  .  (p  (w^,  t) 

V  =  1 

+  V  6*  (w^,.  t)  :  g  (to^.  t)J  .  (2-11) 

Equation  2-11  is  one  way  of  saying  that  the  total  power  flow  is  a  constant  even 
though  it  can  redistribute  itself  amongst  the  several  frequencies.  ABDP  dem¬ 
onstrated  that  for  the  6*  *  f  term  Eq.  2-11  followed  from  the  symmetry  re¬ 
lations.  The  argument  cannot  be  reversed,  however. 

The  physical  argument  that  the  "steady  state"  energy  density  only  de¬ 
pends  on  the  final  state  and  not  on  the  path  by  which  it  was  attained  is  equivalent 
to  requiring  •<8U/ 8t>  .j,  be  an  exact  differential;  i.  e.  ,  OU/8t>,j.  =  d$/dt  , 
or 

n 

d$  =  2  Re  ^  t)-  dn  t)  +  £*  (w^,  t)  •  d  £  (co^,  t) 

V  =  1 

+  V  i*  :  d  g  (w^,  t)]  .  (2-12) 

One  can  define  a  second  potential 

n 

F  =  §  -  2  Re  J  («^,  t)  •  TJl  («^.  t) 

V  =  1 

+  §.*  (w^^  .  t)  •  P  (w^  .  t) 

+  V  6*  (Wy.  t)  :  ft  («^,  t)]  (2-13) 

so  that 

P  (u)„.  t)  =  -  8F/  8  6*  (w  t)  (2-14a) 

t)  =  -  8F/  8^  («^.  t)  (2-14b) 

ft(«^,  t)  =  -  8F/  afrjE*  (u>y.  t)}  . 


(2-14c) 
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From  the  existence  of  a  "time  averaged  free  energy"  and  the  relations  2- 
14a, b,c  one  can  describe  all  of  the  conventional  electro-  and  magneto -optical 
effects  as  well  as  all  of  the  phenomena  by  which  one  obtains  parametric  am¬ 
plification  with  nondissipative  media.  Many  of  the  dispersive  effects  in  mag¬ 
netic  double  resonance  experiments  also  follow.  In  the  subsequent  sections 
we  will  consider  several  forms  for  F  and  the  phenomena  which  they  lead  to. 

For  weak  nonlinearities,  by  the  methods  of  ABDP,  macroscopic  cur¬ 
rent  density  J,  in  Eq.  2-2  and  Eq.  2-3  can  be  partitioned  into  linear  and  non¬ 
linear  parts.  At  the  frequency  co^  one  can  write 

V  xe  (w^)  =  -  i{wyc)  ajc  ((o^)/at  (2-i5) 

V  x3C  =  i{w^^/c)  l{w^)  •  £  (w^) 

+  (4Tr/  c)  ^ 

where  it  has  been  assumed  the  linear  material  is  nonmagnetic:  i.  e.  ,  j£{Wy)=  0 
so  that  =  B  is  given  by  the  Fourier  component 

of  the  nonlinear  part  of  Eq.  2-3.  The  linear  dielectric  constant  at  is  given 
by  €  («^).  An  equivalent  form  of  Eq.  2-15  would  be 

V  X  6  )  =  -i  (w^  /  c)  [X  )  +  4tt  ^  (w^ 

V  (Wy)  =  i  (Wy/ c)  [€  (Wy)  •  6  (Wy)  +  4ir  P  («y)^^ 

-  4irv  •  e  (2-16) 

where 

5e(«-^)  =  ®  (u,y)  -4ir^(Wy)^^. 

The  difference  between  Eqs.  2-15  and  2-16  is  the  difference  in  the  meaning  of 
3C  (w  ).  The  two  different  definitions  of  ?£  (w  )  will  lead  to  different  boundary 
conditions  auid  slightly  different,  although  equivalent,  forms  for  the  energy 
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density  and  power  flow.  Equation  2-15  with  correct  boundary  conditions  will 
prove  most  straightforward  to  use  in  the  approximation  that 
a  constant.  This  is  the  approach  used  by  Bloembergen  and  Pershan  [12]  in 
the  treatment  of  boundary  harmonics.  The  interpretation  of  the  exact  nonlinear 
coupled  amplitude  equations  of  ABDP  will,  be  somewhat  more  direct  for  mag¬ 
netic  nonlinearities  when  Eq.  2-16  is  used. 
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III.  Electric  Dipole  Effects 

We  will  consider  all  effects  that  can  be  derived  from  a  free  energy  of 
the  form  ^  £  6  •  £  &  &  >  ®tc.  ,  as  electrit:  dipole  effects.  The  first  term  is, 

of  course,  the  free  energy  for  a  linear  medium,  but  it  will  prove  illustrative  to 
treat  it  by  the  methods  we  will  use  for  the  more  complicated  phenomena. 
Consider  the  free. energy 

F  =  4  ftj(w)  §i*  (w)  6j(a))  +  §i(a,)  6j*{co)}  (3-1) 

where  it  is  understood  that  repeated  indices,  are  to  be  summed  over.  Since 
there  are  only  nine  independent  products  g^*(w)  6j(w)  there  can  be  only  nine 
independent  tensor  components,  i.  e.  ,  is  Hermetian 

X..(a,)  =  X..*(a))  .  (3-2) 

From  Eq.  2- 14a 

P.(«)  =  -8F/8  6i*(«)  =  )(..(«)  6.(0,) 
and  the  dielectric  tensor  c 

€..(o,)  =  1  +  4ir  X..(o,)  (3-3) 

is  also  Hermetian.  If  the  crystal  is  nonmagnetic,  it  is  invariant  under  time  re¬ 
versal  and  the  tensor  X  must  also  be  invariant  under  time  reversal.  From 
Eq.  A-7  of  the  appendix  one  obtains 

X..(o,)  =  X..*(o,)  (3-4) 

so  that  JC,  and  j,  are  real  symmetric  tensors  as  is  well  known.  Note,  however, 
that  for  magnetic  crystals  the  off  diagonal  elements  are  imaginary  and  one 
obtains  the  well-known  Faraday  rotation.  This  will  be  discussed  at  greater 
length  in  Section  V  when  magnetic  nonlinearities  are  considered.  A  second 
point wdrthnoting,,  however,  is  that  optical  activity  is  not  an  electric  dipole 
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effect,  since  this  requires  imaginary  off  diagonal  elements  of  €  for  a  crystal 
invariant  under  time  reversal  [18],  In  Section  V  it  will  be  shown  that  optical 
activity  follows  from  a  free  energy  of  the  form  X, :  g  3C  . 

The  linear  electro-optic  [19],  or  Pockels,  effect  and  the  largest  of  the 
optical  nonlinearities  can  be  obtained  from  a  free  energy  of  the  form 

+  "2*  "l^  ^i<"3^  ®k*<"l?} 

where 

Since  there  is  only  one  way  to  form  the  product  of  the  ith  component  of 
§*(‘*>2),  the  jth  component  of  £(b>2)>And  the  kth  component  of  §(wj)  ,  there  is  only 
one  quantity  “2’  *^1^  permutation  of  frequencies  and  indices 

must  be  equal,  i.e.  ,  ^2’  *^1^  ~  ^jik^^^Z’  ‘^3’  ‘*^1^’  This  is  not 

yet  the  symmetry  obtained  by  ABDP  which  are  relations  between  the  nonlinear 
polarizability  tensors. 

From  Eq.  2 -14a 

(Pi(«3)  =  ^ijk(‘-3’  "2’  "l>  ^k<“l^ 

=  ’^ijk*<"3'  "2'  “l>  ^i<‘*'3>  ^k*<"l) 

\<“l>  =  ’^ijk*^"3’  "2*  "P  ^i<"3)  • 

The  tensors  of  ABDP,  ,  were  defined  initially  from  the  following  relations, 

o. 


‘^i<"3)  =  ’^ijk<"3’  "2’  "l>  ^k<"l) 

(P.(«2)  =  °’^jik<"2'  "3'  "1^  ^i^"3>  Sk*<"l) 

^«-l)  ”  °’^kij^‘^r  "3’  "2^  6i<“3)  6j*(«2)  • 
From  Eqs.  3-6  and  3-7 


(3-7) 
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‘^2*  =  %k<"3'  "2*  “l)  =  t°V^“2'  "3* 

=  [°^kij^"r  “3*  “2^^*  ' 

For  crystals  invariant  under  time  reversal,  Eq.  A-7  of  the  appendix 

shows 

’^ijk<“3'  “2*  “l)  =  ^ijk*<"3*  “2*  “!>  ’ 
and  Eq.  3-8  is  equivalent  to  the  relations  obtained  by  ABDP.  For  crystals 
which  lack  time-inversion  symmetry,  Eq.  3-8  is  a  generalization  of  those 
relations . 

For  the  particular  case  of  second  harmonic  generation  it  is  simplest 
to  start  from  the  free  energy  rather  than  to  take  the  limits  of  Eqs.  3-6  as 
In  the  limiting  procedure  it  is  difficult  to  keep  track  of  factors 

of  2. 

The  physical  effects  resulting  from  these  types  of  nonlinearities  at 
optical  frequencies  have  been  thoroughly  treated  by  ABDP  and  others  [12]. 

To  obtain  the  low  frequency  linear  Kerr  effect  one  must  take  the  limit 
as  0)^^  0.  The  free  energy,  Eq.  3-5,  must  be  extended  to  include  the  case 

when  =  W2  -  If  is  small  one  should  expect 

Xi.k(«2  +  "r  "2'  "l^  =  V^"2’  "2  -  "l> 

so  that 

F=  -2  Re  X.jj^(w2+Wj,  «2'  “  "l^  ®k^"l^ 

+  §^*(^2^  ®j^"2  ■  "1^  ^k^“l^^  ■ 

If  the  nonlinearity  is  large  enough  so  that  the  amplitudes  of  the  side  bands, 
i.  e.  ,  (>>2  become  comparable  to  the  amplitude  at  free  energy 

must  be  extended  further  to  include  the  other  side  bands;  i.e.  ,  W2  ^ 

<^2  +.  3  <i)  1 ,  •  •  •  ,  etc. 
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IV.  Electric  Quadrupole  Effects 

All  effects  that  are  derivable  from  a  free  energy  of  the  form  X  :  6  V  6  i 
X:  6  g  V  §  ,  etc.  ,  will  be  considered  electric  quadrupole  effects.  This  does 
not  imply  neglecting  electric  dipole  moments  of  the  form  X;  V  J  and  X  :  £  V  £ 
but  rather  that  these  effects  have  their  origin  in  a  quadrupole~type  interaction. 
Restricting  attention,  at  this  time,  to  tensors  X  symmetric  in  the  last  two  in¬ 
dices;  i.e.,  if 

X(u2.  «2’  “j):  £(“2)  V  S  (w  j) 

-  ^jki<"3*  “2*  “1^  6i*(“3)  6j(u)2) 

it  is  assumed 

’^ijki<“3*  ^Z‘  ‘^1)  =  ^ijfk^^^S'  “2'  ‘^1^  • 

That  part  of  X  antisymmetric  in  the  last  two  indices  will  multiply  9j^  ^1^***!)  " 

Ir  Ir  §r\ 

^  i  equal  to  -i(4Tr«j/ c)  «  where  e  is  the  unit 

antisymmetric  tensor  of  the  third  rank.  If  kin  is  a  cyclic  permutation  of  x, 
y,  z,  =  fl;  if  it  is  an  antisymmetric  permutation,  =  _1.  Thus,  the 

part  of  X  antisymmetric  in  the  last  two  indices  can  be  written  as  X :  g  JC  , 

X:  6  €  3Ci  cfc.  Terms  of  this  type  will  be  discussed  in  detail  in  Section  V. 

A  free  energy  of  the  form  X:  6  V  6  i  symmetric  in  the  last  two  indi¬ 
ces,  corresponds  to  a  quadrupole  correction  to  the  linear  dielectric  constant. 
The  third  rank  tensor  X  reverses  sign  on  inversion  of  the  spatial  coordinates: 
i.  e.  ,  X  -►  -X,  y  -*  -y,  z  -♦  -z,  and  thus  vanishes  for  all  crystals  invariant  un¬ 
der  spatial  inversion.  For  the  remaining  crystals  this  correction  has  been 
treated  in  detail  by  Satten  [20]  and  will  not  be  discussed  here.  Terms  of  the 
form  F  =  -  X I  (V  6)(V  £^  are  of  higher  order  and  neglected  here. 
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The  lowest-order  nonlinear  quadrupole  term  has  a  free  energy  of  the 

form 

F  =  -2  Re  “l^  6i*(‘*>3)  6j(w2) 

^fikj^“r  “3’  “2^  6i*(«3) 

+  ^jjfki^“2’  “r  "3^  ®k 

for  0)3  =  cjj  +  ^^2  ’ 

In  addition  to  the  symmetry  of  Eq.  4-1,  one  requires  the  dispersion,  or 
permutation  symmetry  relations  between  the  first  two  indices  of  the  type 

’^ijki<“3'  "2*  “l^  =  ’^jiki<“2*  “3*  “!> 

^iikj<"i*  "3'  "2)  =  ^iikj<"3'  ‘-r  "2> 
etc. 

The  reasoning  behind  Eq.  4-3  is  the  same  as  was  used  to  justify  the  symmetry 
of  the  tensor  in  Eq.  3-5.  The  importance  of  a  free  energy  of  this  type  derives 
from  the  fact  that  for  crystals  invariant  under  space  inversion  it  represents 
the  largest  nonmagnetic  mechanism  for  producing  second  harmonic. 

From  Eqs.  2 -14a,  c  and  4-2  one  obtains 

^i<"3^  =  ^jkf<“3-  "2'  "1)  ®k  ®f<"l^ 

+  "3'  "2>  ®k 

®ki<"3)  =  <jiki<“2'  "r  "3^ 

=  ’^ijki*<"3’  "2’  "1)  ei<"3)  \  ®f*<"l) 

+  ’^jiki*<"2'  "l’  "3^  ®k^i^"3^ 

llkj^***2^  ”  ^iikj*^***l’  ^3*  *^2^  ^l*^***l^  6j(<*>3)  •  (4-4b) 


and 
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^  W^“2*  ‘-1'  ‘"3) 

®ki<"l)  =  ’^ijki*<‘-3*  “2*  “l>  ^i<“3)  ^j*<"2)  ’ 


(4-4c) 


It  is  a  rather  trivial  loss  of  generality  to  restrict  attention  here  to  crystals 
invariant  under  time  reversal.  From  Eq.  A-7  of  the  appendix,  all  the  X's 
are  real  and,  henceforth,  the  stars  can  be  neglected.  The  nonlinear  source 
currents  jjf  ,  to  be  inserted  into  Eqs.  2-15,  are  of  the  form 

,NLS  .  ,  .NLS  .  ^  ,  ,NLS 


=  i  <^3  “2-  "1) 

“  ’^jiki^^Z*  "1*  “3)1  €j<"2)  ®k  M"l> 
‘‘■'“3t*‘iikj<"r  “3*  “2^  ’Siki^^z'  "r  “3)1 


(4-5) 


6/«l)  9k  V"2> 


(4-6a) 


;.(«2)  =  i  <*>2  [^ijkf<"3’  "2’  "l>  ■  ’^fikj<‘^l-  "3*  "2>1 

6i(«3)8j,6/(«l) 

+  '"2tSiki<"2'  "l*  "3>  •’^iikj<"r  "3*  "2^1 


jiki 

i/(«l)  \6i(‘-3) 


(4-6b) 


/|("l)  =  ^  "1  t^iikj<“l'  "3'  "2>  -  ^jkf<"3'  "2*  "l)J 


iikj' 

.ii(<-3)  9j^6j*{u,2) 

"1  t^jiki<"2*  "r  "3>  ■^jki<"3'  "2*  "1^1 

6j*(«2)  ^k  ^i^"3^  •  (4-6c) 


Defining  tensors  X  such  that 
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“2’  "l)  = 

=  i<*>2l^^^(^2’  “P  "3) 

=  icj^  1^^^  (u^,  <1)3,  cj^)  : 

one  can  see  that  the  permutation  symmetry  relations  of  ABDP  (i.  e.  ,  Eqs.  3-8) 

EFF 

do  not  hold  exactly  for  )(  For  example,  if  k-  =  k,  +  k,  +  Ak  ,  from  Eq.4-3 

[Xg^^(u)3,  w^,  Wj)]*  -  "2»  “3^ 

=  i  “2’  “l>  •  "2’  “3)1 

Only  for  the  case  of  exact  phase  matching,  i.  e.  ,  Ak  =  0,  can  one  treat  the  cur¬ 
rent  elements  of  Eqs.  4-6  as  effective  dipole  moments  per  unit  volume.  The 
right-hand  side  of  Eq.  4-8  is  related  to  the  transport  of  energy  by  the  quadrupole 
interactions  as  discussed  in  relation  to  Eq.  2-5.  Since  Jt  is  imaginary,  and 

W3  =  0)^  +  «2 

•  J.(t)^^>=  2Re  ^i«2|x^r^(w2,  Wp  “3)  “  t  ^“3>  "2>  “l^l 

‘"l  <“1>  "3>  "2>  ■  "2’  “1>  1*}] 

£j*(«2)  6.(«3)  . 

Applying  Eqs.  4-3  and  4-8 

<JE{t)- J,(t)^^^>  =  2Re|[i«j  «2’  "P  +  '“2  ^likj^“P  "3»  “2^  i 

■  ^"3’^jlki<"2»"p"3>]  *[<‘‘;^>kV<"l)^j*<"2)®i<“3)]|  > 
and  this  can  be  recognized  as  the  time  average  of  V  •  [  E,(t)  •  »a/8ti. 
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The  coupled  amplitude  equations  developed  by  ABDP  (i.  e.  ,  Eqs.  4-9 
of  ABDP)  can  be  generalized  to  include  quadrupole  nonlinearities  by  replacing 
iw^^  with  •  The  exact  solutions  will  follow  in  the  same 

manner  as  the  dipole  nonlinearity  except  that  the  integration  constant  corre¬ 
sponding  to  the  time  average  Poynting's  vector  will  not  be  simply  <  E  x  ^>,p  , 
but  it  will  have  an  additional  term  corresponding  to  the  transport  of  energy  by 
the  quadrupole  interaction  shown  in  Eq.  4-9. 

The  currents  given  by  Eqs.  4-6  can  be  seen  to  be  invariant  if  one  adds 
to  the  free  energy  (Eq.  4-2)  a  term  that  would  correspond  to  a  surface  energy 
density,  F'  =  F  +  ^  ^ 

^  -  "2-  “l'  ^j<“2>  1 

“’‘ijk<<"3>"2>  "1*  =  “1’  "3*  =  ’‘likj<"l’  "3’  “2*  >  Eq-4-2  reduces 

to  a  surface  energy  density,  and  Eqs.  4-6  all  yield  zero  currents.  It  may  often 

prove  convenient  to  define  new  tensors  by  a  suitable  definition  of  such  that 
F'  is  given  by  Eq.  4-2  with^'  replacing  and 

^<■'"3’  "2’  "l*  ^  ’‘j2ki'<“2-  "l>  "3* 

^likj'^"l’  "3’  “2^  "  °  ■  (4-11) 

Alternatively,  one  could  set  one  of  the  tensors,  for  example,  u^,  «j)=  0, 

with  no  loss  of  generality. 

For  second  harmonic  generation  there  are  only  two  tensors, Jl(«,  u,  2fa>) 
andJ|C(2w,  w,  w) .  With  no  loss  of  generality  one  could  set  X(2w,  u,  w)  =  0  and 
obtain 

vt 

F  =  -  2Re  <d,  2(a)  :  ^(«)  J/w)  ^^(2«) 


(4-12) 
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so  that 

iB(2u>)  =  -  2ia)V.[X(u),  <0,  2(0)  :';§,{a))  gjw)]  .  (4-13) 

For  an  isotropic  material,  the  form  of  «,  2(0)  is  given  by  Table  I.  It  is 

straightforward  to  demonstrate  that  for  a  plane  wave  at  frequency  co  moving 
with  wave  vector  ^(w)  such  that  J^(w)  •  =  0  the  current  given  by  Eq.  4-13 

is  parallel  to  k((o) .  Neglecting  surface  harmonics  [12]  this  current  cannot 
generate  a  second  harmonic. 

If  the  field 6((o)  is  not  a  simple  plane  wave,  but  two  plane  waves  as 
shown  in  Fig.  II,  it  is  possible  to  match  phase  velocities  in  the  z-direction,  if 
the  isotropic  material  exhibits  anomalous  dispersion.  The  sum  of  ^(w)  and 
^'((o)  has  components 

6x(w)  =  ^  M  ^11  '  1  cosk^x-i[6||-§||*]  cos  4  sink^x^  x  exp  (-  ik^  z  ) 

6y(«)  =  '  ]  cosk^  X  -  i[6j^  -  8^^']  exp(-  ik^  z) 

=  |[  6  -  6  ,,  ]  sin<^  cos  X  +  i  [  6  II  '  +  6  „  ] 

sin(^sink^x^  exp(-ik^z)  (4-14) 

The  quadrupole  moments  and  the  currents  can  be  obtained  by  direct  substitution 

of  Eqs.4-14  into  Eq.  4-13  where  ^is  given  in  Table  I.  The  only  terms  that 

will  couple  to  a  plane  wave  at  2(u  propagating  in  the  z-direction  are  the  parts  of 

fi  and  Q  that  do  not  have  an  x- spatial  dependence.  The  term  fi„_(2(i)) 

"^xz  "^yz  ^  xz 

vanishes  identically,  however, 

=  ^6,  6  ®  .1 

^^NLS(2«)  =  -  4«k^  -  g„  8^']  sin<^exp(-i2k^z)  .  (4-15) 
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Observe  that  if  there  is  symmetry  between  £  and  such  that  g||  '  =  6j|  and 
gj_  =  5j_'  ,  there  is  no  second  harmonic  generation.  This  effect  can  be  large 

2  1/2 

only  for  large  anomalous  dispersion  since  sin  4>  =  [1  >  (n(2u>)/n(a>))  ]  '  . 

Second  harmonic  generation  by  means  of  a  quadrupole  nonlinearity  will 
thus  be  a  very  weak  effect  in  all  isotropic  media. 

Anisotropic  media,  on  the  other  hand,  can  have  observable  second  har¬ 
monic  generation  by  a  quadrupole  nonlinearity.  Consider  the  tensor  cj,  2u) 

for  calcite  given  in  Table  II.  There  are  several  differences  between  this  and 
the  fourth  rank  elastic  tensors.  For  example,  K  (cu,  w,  2u>)  /  X  ((d,  w,  2u), 
X^y^(w,  cj,  ’  ****  because  there  are  no  operations  that  will 

transform  z  into  x  or  y  and  leave  the  crystal  invariant.  The  identity 

X  (hi,  cj.  2(1))  =  X  (w,  (>),  2<t))  follows  from  the  operations  of  a  trigonal 
xxyy  yyxx 

axis  in  calcite. 

One  can  match  the  phase  velocity  of  an  ordinary  ray  at  the  fundamental 

frequency  to  an  extraordinary  ray  at  the  second  harmonic  [  3,  4]  as  shown 

in  Fig.  III.  As  discussed  in  ABOP  (£q.  4-8),  harmonic  generation  is  due  to  the 
[NLS, 


component  of  ^  (2(>))  parallel  to ^<6., 


'useful 


=  -(4fa)^  c  ^(w)/ c)  6^(w) 


f^^2^x  Vy'z'x'^"’ 

^®2^z  ’^y'y'z'z'^"'  ’  (4-16) 

Although  the  linear  optical  propestiea  of  a  uniaxial  crystal  are  constant  for  all 
rays  on  a  cone  forming  a  given  angle  0  with  the  crystal  axis,  this  is  not  true  for 
the Tiohlinear  properties.  The  most  general  form  of  the  tensor  components  in 
Eq.  4-16  can  be  obtained  from  Table  II  by  a  rotation  through  an  angle  (|>  about 
the  z  axis  (i.  e.  ,  crystal  axis)  so  that  the  new  y'  axis  is  parallel  to  6  (<■))  and 
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then  a  rotation  by  0  about  g,(cj)  so  that  z'  is  the  direction  of  h(w). 

.s  rfNLS,,  X  9  2  1/2,  ,  -1  -2,  . 

®2  *  i  (2«)  =  -2a)  c  '  (oj)  c  6  (u) 

X  2  +  ^13)  O’  +  3  “  2^ 

+  2  ^  sin  3<|>  cos  (20  -  a)^  ‘  (4-17) 

The  fact  that  Eq.  4-17  is  non -zero  is  due  to  the  crystal  anisotropy  since  for 

isotropic  crystals  the  primary  wave  is  a  transverse  wave  and  a  =  0.  From 

Table  I  isotropic  crystals  have  ^22=  ^^nd  4=  0- 

The  experiments  of  Terhune  et  al.  [5]  detected  quadrupole  second  har^ 

monic  in  calcite  for  one  particular  geometry.  They  calculate  for  their  orien- 

- 1 8 

tation  of  crystal,  the  bracketed  term  in  Eq.  4-17  was  of  the  order  of  10 

esu  units.  This  sets  a  lower  limit  on  the  three  constants  X  ^  2>^]^  3'  4 

since  <j>  might  have  been  set  near  0,  ■jr/  3,  2ir/  3,  etc.  ;  the  real  X 's  may  be 

larger  by  approximately  (1/ sin  a)^  20.  It  would  be  interesting  to  see  if  an 

angular  dependence  corresponding  to  the  last  term  is  observable.  This 

would  give  a  direct  measure  of  X,  ^  independent  of  X,  _  and  X,  -  . 

1,4  ^  1,2  lj3 

Although  in  isotropic  material  there  is  no  second  harmonic  generation 
per  unit  volume,  due  to  the  quadrupole  nonlinearity  there  can  be  generation 
at  the  surface  [12].  The  quadrupole  term  in  the  current  density,  i.e.  ,  Eq. 

4-5  changes  the  boundary  condition  that  the  tangential  component  of  ^  is  con¬ 
tinuous  across  the' surface  of  a  nonlinear  dielectric.  One  can  show  from 
Eq.  2-15  that  if  2  ^  vector,  normal  to  the  surface  between  medium.  1 

and  medium  2,  directed  from  1  into  2;  and  if  is  a  unit  vectdr  parallel  to  the 
surface,  such  that  t^i  7’  ^  =  0,  the  boundary  condition  on  3C(a))  can  be  written 

-5Ci(«)).  -(4iri<.,/c)  [e2(«)  -)ei(w)]:%^2(^l,2*^)‘ 


(4-18) 
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If  the  fundamental  is  polarized  perpendicular  to  the  plane  formed  by  the  normal 
to  the  interface  and  the  direction  of  wave  propagation,  Eq.  4-18  reduces  to  the 
usual  condition  that  the  tangential  component  of  ^  is  continuous.  In  this  case, 
the  quadrupole  surface  harmonic  is  obtained  by  substituting  for  of  Eq. 

4-12  of  reference  12,  the  quantity  (iw)  ^  as  given  by  Eq.  4-13  of  this 

paper.  The  tensor  is  given  in  Table  I. 

For  other  polarizations  of  the  fundamental,  the  right-hand  side  of  Eq. 
4-18  does  not  reduce  to  zero  and  to  calculate  the  boundary  harmonics  it  must 
be  used  instead  of  Eq.  4-11  of  reference  12. 
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V.  Magnetic  Dipole  Effects 

All  effects  which  can  be  derived  from  a  free  energy  proportional  to 
one  or  more  powers  of  the  magnetic  field  will  be  considered  magnetic  dipole 
effects.  Higher  magnetic  multipoles  are  explicitly  neglected  since  free  en¬ 
ergies  proportional  to  gradients  of  the  magnetic  field  are  not  treated.  The 
simplest  magnetic  free  energy  describes  the  linear  phenomena  of  optical 
activity  [18], 

F  =  -  {x..(co)  §.*(0,)  3Cj(«)  +  )(...*  {«)  §.(w)  3Cj*  (u))}  .  (5-1) 


Assuming  a  crystal  invariant  under  time  reversal,  Eq.  A-7  of  the  appendix 
requires  )(-(<»))  =  -  (u);  i.e.  ,  X  is  pure  imaginary.  The  forms  of  X  for 

quartz  and  NaClO-j,  both  optically  active,  are  given  in  Table  III.  The  die¬ 
lectric  currents  are  obtained  from  ^  (a)  =  80(u))/  9t  +  c  V  x  and  one 

can  write  the  part  of  the  dielectric  constant  due  to  Eq.  5-1  in  the  form 

fin 


A  €i^  =  i  4ir  (c  ky  u)  (lX..|  +  |X^|)  e' 


(5-2) 


fin  . 


where  k  is  the  propagation  vector  for  the  wave  and  c  is  the  antisymmetric 
third  rank  tensor  introduced  in  Section  IV.  The  second  rank  pseudotensor  X 
willi  vanish  if  the  crystal  has  inversion  symmetry  since  the  second  index 
transforms  like  JC  and  keeps  its  sign  under  inversion.  The  somewhat  weaker 
condition  of  the  presence  of  a  mirror  plane,  however,  can  make  individual 
terms  in  X  vanish.  For  example,  if  the  crystal  is  invariant  under  reflection 
in  the  x-y  plane,  the  only  nonvanishing  elements  of  X  are  X^^, 

X^y,  This  has  the  physical  significance  that  if  light  is  propagating  parallel 
to  a  crystal  mirror  plane,  the  crystal  must  have  the  same  effect  on  right  and 
left  circular  polarization  and  there  can  be  no  optical  activity.  Note  that  the 
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sign  of  the  effective  dielectric  constant  in  Eq.  5-2  depends  on  the  sign  of  k^. 
This  has  the  important  consequence  that  light  going  forward  and  backward 
through  an  optically  active  material  will  emerge  with  no  net  rotation  of  the 
plane  of  polarization.  This  should  be  distinguished  from  the  Faraday  effect 
which  doubles  the  rotation  on  passing  the  same  crystal  backward  and  forward. 

The  form  of  free  energy  leading  to  the  Faraday  effect  represents  the 
simplest  type  of  magnetic  nonlinearity, 

F=  -2  Re  [  ‘*>2*  6j(w2)  JCj^(w  j) 

^jki^‘^2’  “l’  “3^  6j(«2^  ^i*^“3^ 

+  ^kij^“r  “3’  “2^  ®k^"P  §i*(«3)  3Cj(w2)]  (5-3) 

where  co^  =  Considering,  for  the  moment,  only  crystals  invariant 

under  time  reversal,  Eq.  A-7  of  the  appendix  requires  these  X's  to  be  pure 
imaginary.  From  arguments  used  in  both  Sections  III  and  IV,  there  is  a 
permutation  symmetry  between  the  first  two  indices 

Xijk(«3.  «2»  "i)  =  ^jik^"2’  "3’  "l^ 
etc. 

If  Wj  <<  W31  orie  must  also  include  the  free  energy  for  co^'  =  “  Wj  in 

the  same  manner  as  was  discussed  in  Section  III  for  the  linear  electro-optic 
effect.  Let  « ^  correspond  to  a  low  frequency  (i.  e.  ,  w  j/  2ir  <  10^  cps)  while 
a^nd  <ti^  are  optical  frequencies.  From  the  normal  dispersion  of  magnetic 
phenomena  [21],  it  follows  that  Mf^2.  i.“l’  “2’  “l^  much  larger  thaui  the 
other  four  tensors  and  one  can  simplify  the  free  energy  by  taking 

3Cijk(«2+  "l»  "2>  ‘^1)  ^ijk<"2»  “2  ■  "1’  "1)  =  V 

F  =  -2  Re  [6i*(«2  +  "i>  €j(«2)  3ICk(«  1) 


(5-5) 
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One  obtains 


(P.(a)2  +  wj)  =  g.(«2)3Cj,(u>i)  (5-6a) 

(p.(a;2  -  6j(«2)3ej^*(a>l)  (5-6b) 

^  gj*(w2  -Wj)].  {5-6c) 

In  the  limit  that  -*•  0,  Eqs.  5-6a  and  5-6b  will  describe  the  dc.  Faraday 
effect.  One  should  not  interpret  the  symmetry  of  Eq.  5-4  to  mean 
Xijk(w.  w,  0)  =  since  both  terms  in  Eq.  5-4  multiply  («)  6j(w) 

in  the  limit  that  -*  0  in  Eq.  5-3. 

For  an  isotropic  material,  invariant  under  time  reversal,  one  can  show 
the  only  non-zero  elements  of  X  are 


X  =X  =X  =-X  =-X  =-X  =i|X,  (5-7) 

xyz  yzx  zxy  xzy  zyx  yxz  '1,^3'  ' 


and  one  can  write  Eqs.  5-69.,  b  as 

?(«2)  =  i  I’^i,2,3l6<“2)  + 


(5-8) 


The  change  in  the  effective  dielectric  constant  thus  has  the  usual  form 

A*..  =  i4Tr  |Xi^^3l€'j‘^hj^(t)  (5-9) 

where  hj^(t)  is  the  real  value  of  the  kth  component  of  the  magnetic  field  at  Uj^. 
The  form  of  Eq.  5-9  is  identical  to  the  form  of  Eq.  5-2  except  that  the  latter 
changes  sign  with  reversal  of  the  direction  of  propagation. 

If  one  initially  has  a  circularly  polarized  optical  field  propagating  in 
the  z -direction,  Eq.  5 -6c  shows  there  will  be  a  z -component  of  magnetization 
at  zero  frequency  [22]. 


»z«»  =  ±2  l*l,2,3ll«l^  • 

In  concentrated  neodymium  ethylsulphate 


known  from  Faraday 


(5-10) 
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*  9  o 

measurements  to  be  approximately  10  cgs  units  at  4,  2  K  [23],  A  one-mega- 

r  -4 

watt  laser  pulse  will  thus  produce  a  total  flux  §  =  \  4ttM  •  dA  =  10  cgs  units. 

If  the  optical  pulse  is  lO"^  seconds  long  and  is  detected  by  a  1000  turn  coil, 
the  pick-up  voltage  should  be  1  millivolt.  This  is  independent  of  the  cross- 
section  area  of  the  laser  pulse  so  long  as  the  volume  over  which  M  is  produced 
is  large  enough  that  the  flux  does  not  close  on  itself  completely  within  the  pick¬ 
up  coil. 

The  nonlinearity  described  by  Eq.  5-3  is  of  no  importance,  for  media 
invariant  under  time  reversal,  when  becomes  an  optical  frequency.  Mate¬ 
rials  which  lack  time -reversal  symmetry,  either  because  they  are  subjected 
to  an  external  dc  .  magnetic  field  or  because  they  are  ferromagnetic,  can  have 
detectable  nonlinear  effects  even  when  co  ^  is  an  optical  frequency.  Consider  a 
cubic  crystal,  class  that  lacks  time -inversion  symmetry  because  it  has  a 
dc  .  magnetization  in  the  (001)  direction.  The  point  group  operations  under  which 
this  crystal  is  invariant  are  shown  in  Fig.  IV.  The  form  of  is 

given  in  Table  IV.  If  the  crystal  were  invariant  under  time -reversal  symmetry, 
the  real  parts  of  the  tensor  given  by  Table  IV  would  vanish  and  |  |  would 

equal  |  as  for  the  dc,  Faraday  effect  described  by  Eq.  5-7.  The  real 

terms  can  be  seen  to  come  from  the  dc .  magnetization  in  the  (001)  direction 
by  considering  a  cubic  crystal,  invariant  under  time  reversal,  but  subjected 
to  a  dc,  magnetic  field  in  the  (001)  direction.  The  free  energy  will  have  terms 
of  the  form 

F  =  -2  Re  ^iji^(<*>3.  «2’  "l^  6.*(«3)  Gj(w2)5Cj^(w j) 

-2  Re  «2*  "r  6i*(«3)  6j(«2^^k^" * 


(5-11) 


FIG.  4.  THE  POINT  GROUP  SYMMETRIES  OF  A  CUBIC  CRYSTAL,  Oj^  , 
WHICH  HAS  A  DC  MAGNETIZATION  ALONG  THE  (001)  AXIS. 
FIGURE  4a  SHOWS  THE  PURE  SPATIAL  OPERATIONS  UNDER 
WHICH  THE  CRYSTAL  IS  INVARIENT.  FIGURE  4b  SHOWS  THOSE 
SPATIAL  OPERATIONS  WHICH  WHEN  COUPLED  WITH  TIME 
REVERSAL  WILL  LEAVE  THE  CRYSTAL  INVARIENT. 


TR393 


-26- 


Since  the  crystal  is  invariant  under  time  reversal,  Eq.  A-7  of  the  appendix 

requires  ‘*'1)  pure  imaginary  and  of  the  form  of  Eq.  5-7. 

The  components  ^2’  *^1’  real  and  have  the  same  form  as  the 

tensor  given  by  Table  I  except  that  due  to  dispersion  =  ^zyyz  “^66^ 

X  =  X  =  X/  Similar  relations  hold  for  the  other  terms.  The  free 

zyzy  yzyz  o,  o 

energy  in  Eq.  5-11  can  be  put  in  the  form  of  Eq.  5-3  by  contracting  on  the 

index  i.  The  third  rank  tensor  will  thus  have  the  form  in  Table  V.  If  |  X 
.1  .  r(3) 


4.1 


I  ^6  3^  ”  ”  ^  ^4  1'  ^  ”^12  3^®'^^®®  ^  have  identical  forms.  The  Voigt 

effect  [24],  microwave  modulation  by  the  Faraday  effect  [25,  26],  and  the  non¬ 
linear  optical  effects  correspond  to  the  tensors  in  Tables  IV  and  V  for  =  0, 
near  a  microwave  resonance  frequency,  and  an  optical  frequency,  re¬ 
spectively. 

Consider  the  case  where  3C.(0) is  a  dc  field  along  the  z-axis,  is 

a  circularly  polarized  microwave  field  in  the  x-y  plane,  is  near  a  magnetic 
resonance,  and  light  is  propagating  in  the  x-direction.  Analogous  to  Eq.  5-9 
one  obtains 

A  =  4ri  I  |:h^(t)  t  4,  |  je^(O)  .  (5-lZ) 

Bloembergen,  Pershan,  and  Wilcox  [25]  have  shown  that  for  light  propagating 
parallel  to  the  magnetization  of  a  sample  there  is  a  Faraday  rotation  propor¬ 
tional  to  the  instantaneous  magnetization,  even  if  that  magnetization  is  rotat¬ 
ing  at  a  microwave  frequency.  Mathematically,  a  magnetization  in  the  x- 
direction  will  produce  a  change  in  A  =  i  K  where  K  is  a  function  of  the 
material  and  the  wavelength  of  the  light.  Neglecting  damping,  the  Bloch 
equations  can  be  solved  for 

=  y  yno(«o^  "  <"o  ^x  " '  "i  V 
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and 


A  =  K  V  ‘  +  “l  V 


Comparing  Eqs.  5-12  and  5-13 

.(3)  ,  „  2 


4tr  |x\^y  =  K  7n^3C^(0)  (a,/  - 


2,-1 


4ir  =  K  Y  uij 


,  2  2.-1 

“i'“o  ■  "l  ' 


(5-13) 

(5 -14a) 
(5-14b) 


where  7>l^  is  the  steady  state  dc  .  magnetization,  is  the  microwave  reso¬ 
nance  frequency  Y3C  ^(0),  and  is  the  static  susceptibility  7n^/3C^(0), 

/  O  \ 

In  the  limit  cj,  -*  0,  Eq.  5-14a  approaches  4it  |X,  ,  =  K  X , 

*  qc  ‘ 

1. 3  X  lO"®  cgs  units  in  concentrated  neodymium  ethylsulphate  at  4.  2°K,  this 

being  known  from  Faraday  effect  measurements.  For  cjj«5  co^,  I  X ^^2  3!*^ 

I  (4)  I  *  * 

I  X^  ^1  and  the  discussions  of  Bloembergen  et  al,  on  modulation  of  light  follow. 

(3)  *2 

The  nonlinear  optical  problem  is  obtained  for  >>  w^,  4it  |  X^  ^  3 1  «  ”  <%»  0, 

4ir  »  yK  X^c  ruby  laser  line  y/<*>j  *  10  cgs  units 

r  (4)  f  -14  4 

and  4ir  |X^  '^1%  9  x  10  cgs  units.  In  an  external  dc .  field  of  10  gauss, 

the  nonlinear  polarization  at  will  be  given  by 

EFF 

<Py(«3)  =  («3'  "2*  "l^  6z<"2>3Cy<"l) 

EFF  -10 

where  X  ~  0.  7  x  10  cgs  units  in  concentrated  neodymium  ethylsulphate 
at  4.  2°K.  This  should  be  compared  with  X'  3  x  10”^^  in  KDP  [5]  for  the  elec¬ 
tric  dipole  nonlinearity  at  room  temperature.  Neodymium  ethylsulphate  is  not 
an  isotropic  crystal, but  the  essential  features  of  the  nonlinear  effect  and  the 
estimate  of  its  size  are  not  affected  by  this. 

The  coupled  amplitude  equations  of  ABDP  can  be  obtained  by  replacing 


,NLS 


dNUS 


NLS, 


NLS, 


(w  )  with  g  (w  )  =  i  «  ^  (w  )  +  c  V  x*g,  {u  ),  ^ 


sNLS 


and 


% 


NLS 


following  directly  from  Eqs.  2-14  and  the  assumed  form  of  the  free 
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energy.  In  this  case,  the  boundary  conditions  [12]  at  the  surface  of  the  non- 
linear  dielectric  are  changed  to  the  tangential  component  of  3C  -  4iTTn 
is  continuous  rather  than  just  the  tangential  component  of3C. 

An  alternative,  but  completely  equivalent,  procedure  would  be  to  re¬ 
define  the  quantities  in  Eq.  2-15  so  that  =  3C  (w^)  - 

Neglecting  the  quadrupole  terms,  Eq.  2-15  becomes 

V  x;^(cj^)  =  i  (w^/  c)  €  •  6(w^)  +  4Tr  (i  q^/  c)  • 

(5-15) 

In  this  form  the  tangential  component  of3C'(u^)  is  continuous. 

For  the  exact  solutions  to  the  nonlinear  coupled  amplitude  equations, 
one  of  the  integration  constants  corresponds  to  the  power  flow  being  constant. 
The  proper  form  of  this  term  will  automatically  follow  from  the  equations; 
however,  one  should  note  that  this  constant  will  correspond  to  2Re^l^(w^)x 

V 

using C(C'  rather  thanSC.  This  has  been  pointed  out  in  Section  11  and  is 
analogous  to  the  considerations  in  Section  IV  for  the  power  flow  by  means  of 
a  quadrupole  nonlinearity. 
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VI.  Conclusion 

The  main  purpose  of  this  paper  has  been  to  show  that  for  nondissipa- 
tive  media,  there  exists  a  function  F,  the  time  average  free  energy,  from 
which  one  can  derive  all  the  constitutive  equations  involved  in  the  electro¬ 
magnetic  theory  of  macroscopic  media.  This  includes  linear  as  well  as  non¬ 
linear  relations.  Linear  and  quadratic  electro-optic  effects,  Faraday  and 
magnetic  Kerr  effects,  optical  activity,  as  well  as  the  new  nonlinear  optical 
phenomena  of  harmonic  generation,  mixing,  electric  rectification  and  mag¬ 
netic  rectification  have  all  been  derived  phenomenologically  from  several  of 
the  simplest  possible  forms  for  F.  In  addition,  several  of  the  phenomena  in 
different  frequency  ranges  have  been  shown  to  be  related.  It  has  thus  been 
possible  to  predict  the  order  of  magnitude  of  several,  as  yet,  unobserved 
effects. 

In  principle,  there  is  no  reason  why  these  methods  cannot  also  be 
used  to  define  free  energies  that  are  functions  of  acoustic  fields  as  well  as 
the  products  of  acoustic  fields  and  electromagnetic  fields.  In  this  way  one 
can  also  obtain  the  linear  and  nonlinear  electro-elastic  and  magneto -elastic 
effects.  One  shovild  also  be  able  to  obtain  information  on  dissipative 
effects  by  introducing  several  general  types  of  phenomenological  loss  terms. 
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Appendix:  Time-Reversal  Transformation 
Consider  a  vector  or  pseudovector  quantity  A(t)  and  its  Fourier  trans- 
formA(cj) 


A(t)  =  ^  exp  (iwt)  dear  (A-1) 

-00 


A(w)  =  ^  j0.(t)  exp  (ricot)  dt/2Tr  . 

-oo 


If  A(t)  is  real,  A{u)  =  A*(-a)).  Under  time  reversal,  t  goes  into  -t  and  w  goes 
into  -«  so  that  if  Tj^  is  the  time-reversal  operator 


Tj^  =  J  [Tj^  A(t)]  exp  (-iwt)  dt/2ir  . 


(A-2) 


-00 


If  Tj^  A(t)  =  A(-t)  as  it  does  for  E(t) 


J  A(-t)  exp  (-i«t)  dt/2iT  =  ^(u»)  . 
-00 


(A-3) 


If  Tj^  A(t)  =  -  A(,^t)  as  it  does  for  H(t) 


Tr  ^(w)  =  -  ^  exp  dt/  2ir  =  -^(w)  . 

-oo 


(A-4) 


Consider  a  real  scalar  quantity  $ 

«  =  2  Re  {X  •  •  } 

(A-5) 

where  each  of  the  quantities  A,  B,  C  is  either  a  vector  or  a  pseudovector. 
Under  time  reversal  A  will  transform  like  T„  4l(«  ).  =  t .  where  t  * 

Xv.cIIa  a 
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is  either  I- 1  if  A  transforma  like  ^  or  -  1  if  A  transforms  like 
transformations  hold  for  ^  ^  etc.  Then,  under  time  reversal 

fR  *■  *•  (tr  X(-..  “b.  “c-  •  >ljk. . . 

a* ( c* ^ ^ 

where  n  is  the  nximber  of  quantities  JSt*  transform 

If  ft  is  to  b«  a  real  scalar,  X  must  transform  as 

’’r  i<“a»  “c»  •••)=<■  1  “c‘  •  •  *  ^ 


Similar 

(A-6) 
like  H. 

(A-7) 


where  the  tensor  X  has  n  indices  that  transform  like  ^  on  time  reversal- 
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Table  I 

The  form  of  the  tensor  \(u),  to,  Zw )  to  be  used 
in  Eqs.  4-12  and  4-13  for  isotropic  materials 


XX 

yy 

zz 

yz 

zy 

zx 

XZ 

xy 

yx 

XX 

Xl,l 

Xl,2 

Xl,2 

yy 

Xl,2 

Xl,l 

Xl,2 

zz 

Xl,2 

Xl,2 

Xl,l 

yz 

Xb,  6 

Xb,  b 

zy 

Xb,  b 

Xb,  b 

zx 

Xb,  b 

Xb,  b 

xz 

Xb,  b 

Xb,  b 

xy 

Xb,  b 

Xb,  b 

yx 

Xb,b 

Xb,b 

where 

^6.  6  ^ 

txi, 

.1  ■ 

Xi  2^ 

f 
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Table  II 

In  the  form  of  the  tensor  xiu,  u,  2(J )  to  be 
used  in  Eqs,  4-12  and  4-13  for  calcite 


ijV 

XX 

yy 

zz 

yz 

XX 

Xl,l 

Xl,2 

Xl,3 

Xl,4 

yy 

Xl,2 

Xl,  1 

Xl,3 

■Xl,4 

zz 

X3,  1 

X3, 1 

X3,3 

yz 

X4,  1 

■  X4, 1 

X4,  4 

zy 

X4,  1 

‘X4, 1 

X4,  4 

xz 

zx 

xy 

yx 

where 

Xb,  6 

-r  Ixi  1  -  X, 

z! 

zy  xz  zx  xy  yz 

Xl,  4 

'^1,4 

4 

^4,  4 


X4,  4 

X4,  4 

X4,  1 

X4,  1 

X4,  4 

X4,  4 

X4,  1 

X4,  1 

Xl,4 

Xl,4 

X6,  6 

X6,  6 

Xl,4 

Xl,  4 

X6,  6 

X6,  6 
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Table  III 

The  form  of  the  pseudotensor  x  to  be  used  in  Eq.  5-1  to  describe 

'.iMr 

optical  activity  in  quartz  and  NaClO^.  For  NaClO^  1  ”  ^3  3 

The  symbol  |  Xj  j  |  should  be  taken  to  mean  the  imaginary  part  of 
Xj  ^  Sind  is  thus  a  real  ntimber  than  can  be  either  positive  or  negative. 


X 

y 

z 


ilXj^ll 


1 1X3^  3  I 
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Table  IV 

The  form  of  the  pseudotensor  x(^3>  ^2’  used  in 

Eq.  5-3  for  a  cubic  crystal  that  lacks  time-reversal  sym- 
metry  because  of  a  dc  magnetization  in  the  (001)  direction 


Xijk<“v  “2>  “l' 


ij/k 

X 

y 

z 

XX 

1  ’‘l ,  3  1 

yy 

Ixi.jl 

zz 

1  Xj,  3  1 

yz 

i  1  X4_  1  1 

1  X4^  2  ! 

zy 

i  1  1  1 

1  X4_  2  1 

zx 

1  x4^2  1 

-i  1  X'4,  1  1 

xz 

-*  1  X4, 1  1 

xy 

*  1  X<,^3 

yx 

1  Xe,  3 

Note  that  the  s  I  are  the  real  and  imaginary  parts  of  X.  .  and  can  be 


either  positive  or  negative. 
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Table  V 

The  form  of  the  pseudotensor  that  results  from  contracting  the  free  energy 
of  Eq.  5-11  on  the  index  f.  The  superscripts  3  and  4  refer  to  terms  that 
originate  from  the  third  rank  tensor  and  fourth  rank  tensor,  respectively 


ij/k 

X 

y 

z 

XX 

1  x[%  |3e  (0) 

yy 

1  xS^>2  |3C  (0) 

zz 

1  Xi‘^1  13^  0) 

yz 

‘  1  >‘1!  3  1 

■  1X6%! 

zy 

■*  1  4^2,  3  1 

zx 

‘  1  Xj^2_  3  1 

xz 

xy 

‘  1  3 

yx 

2,  3 

*Note  that  (  and  i  real  and  imaginary  parts  of 

j  and  ^1^2  2  z’op^cf^'v^cly  and  can  be  either  positive  or  negative. 
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